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Abstract. Using a gap equation for the pion mass a nonperturbative method is 
ij_) | given for solving the chiral quark-meson model in the chiral limit at the lowest 

_C ■ order in the fermion contributions encountered in a large N/ approximation. 

The location of the tricritical point is analytically determined. A mean field 
potential is constructed from which critical exponents can be obtained. 
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1. Introduction 

Due to the possibility of exploring parts of the phase diagram of strongly inter- 
acting matter in current (RHIC) or forthcoming (LHC, GSI) experiments, both 
numerical (lattice) and analytical investigations have recently received much atten- 
tion. Theoretically, the phase boundary of QCD depends on the number of quark 
flavors taken into account and the masses of the quarks [ 1 . At vanishing chemical 
potential QCD with two massless quarks undergoes as a function of temperature a 
second order transition from a hadron to a quark-gluon dominated phase. At zero 
temperature and large values of the chemical potential the transition from hadronic 
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to more exotic phases is of first order. With increasing temperature, the line of the 
first order transitions ends in a tricritical point (TCP). For non- vanishing masses 
of the u and d quarks TCP becomes a critical end point (CEP) which separates the 
first order line from a crossover line defined by the locations of finite maxima of 
some susceptibilities. The decreasing value of the strange quark has the same effect 
as increasing the value of the chemical potential, namely it turns the transition into 
a first order one. Therefore, one can expect that for a physical value of the strange 
quark mass the CEP in the \i — T plane will be closer to the temperature axis. 

To locate TCP/CEP nonperturbative methods are required. First principle 
lattice investigations using MC methods based on importance sampling are made 
useless by the fact that the fermion determinant in the partition function is complex 
at non-vanishing baryonic chemical potential. To avoid this problem, recently new 
techniques were developed (for a review see e.g. [ 2 ) to extract, in the /i/T < 1 
region, information for finite /i from simulations at zero chemical potential (multi- 
parameter reweighting of observables, Taylor expansion of the fermion determinant) 
or at imaginary chemical potential, followed by analytical continuation with or with- 
out Taylor expansion. Using the multiparameter reweighting method the CEP was 
obtained for rif = 2 + 1 dynamical staggered quarks with physical masses in [ E] . 
The location of CEP was estimated in the same range in [|1| using an extrapolation 
based on simulation performed with 3 degenerate quarks. 

Analytically, the phase structure of QCD and the location of TCP/CEP can be 
obtained using the Schwingcr-Dyson approach [[5] or the strong coupling expansion 
of the lattice action [|Sj. Alternatively, one can investigate effective models of QCD 
[{7\. Due to its low energy nature, an effective model might prove not satisfactory 
in giving the exact location of the endpoint and the shape of the phase boundary, 
but one can still expect to obtain some insight into the physics near TCP/CEP 
such as noncquilibrium dynamics. It can tell us what is the soft mode at the CEP, 
and by predicting specific signatures [|H] it may help finding its location. Traces of 
a possible overlap between the critical regions of TCP and CEP manifesting itself 
in the \x — T plane as discussed in [ |§] could be in principle verified experimentally 
since it would appear as a gradual change in the value of the critical exponents as 
the systems pass closer and closer to CEP. 

The outline of this paper is as follows. In sec. |21we introduce an effective quark- 
meson model and present our method for its solution. In sec. OHwe summarise the 
main results obtained for the phase diagram and investigate the relation between 
the effective potential of the theory and the conventional treatment of the TCP in 
terms of a mean field Landau-type potential. We conclude in sec. 0J 

2. The model and the large Nf method of its solution 

In the broken symmetry phase characterised by the homogeneous vacuum expecta- 
tion value $ the usual shift in the radial direction a — > yNQ+a and the requirement 
of a finite constituent quark mass m q = g$ as N — v oo defines the following effective 
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SUl(Nj) x SUn(Nf) symmetric Lagrangian for N mesonic and Nj fcrmionic fields 
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counterterm Lagrangian. For fermions the chemical potential is introduced by the 
replacement: dt — > dt~ ip- One can see that the fermions start contributing at level 
0(l/\fN) in the large N expansion, that is at an intermediate level between LO 
and NLO of the mesonic sector. Since we want to calculate fermionic corrections 
to the LO which has O(N) symmetry we disregard the second independent quartic 
coupling of the SU(Nf) x SU(Nf) linear a model which is present for Nf > 2 and 
is proportional with the product of two totally symmetric structure constants d. 

Our method consists in taking into account perturbatively Nf = y/N + 1 flavors 
of N c coloured quarks to 0(l/y/N) at one loop order both in the EoS 



V' ff (<f>)/VN =: $iJ($ 2 ) = $ 



-4 + -T B (M) - - r T F (m q ) 



0, (2) 



and in the gap equation iG~ l (p 2 = M 2 ) = which resums the superdaisy diagrams 
made of pions with one possible one loop fermion insertion. T B i F \ is the bosonic 
(fermionic) tadpole contribution with mass M(m q ). The inverse pion propagator 
reads as 

iG-\p) =p 2 - m 2 - ^ 2 - ~T B (M) - B 14> (p). (3) 

The contribution of the fermionic one- loop bubble ("fish" diagram) B\^ can be 
related to the fermionic tadpole: 



Bii,(p) = - 



4g 2 N c 



T F (m q ) + 



2g 2 N c 



p 2 I F (p,m q ), 



(1) 



where I F (p, m q ) is the scalar bubble integral with external momentum p evaluated 

with mass m q and with Fermi-Dirac distribution, f F (u>) = l/(exp(/3(ui =F //)) + 1). 

Substituting relation (J3J into (0 and making use of (J2J one finds that for p 2 = the 

inverse propagator vanishes, that is Goldstone's theorem is fulfilled at the minimum 

of the effective potential i.e. M = 0. Another way for seeing this is to use the 

definition of iJ($ 2 ) given in (J5J in the expression of the gap equation. Then one 

has 

2g 2 N c 



H(<P Z ) = M z 1 - 
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which shows that in the spontaneously symmetry broken regime defined by iJ($ 2 ) = 
0, $ ^ the physical pion mass is always zero. 

3. Analytical determination of the TCP 

One can choose to study the phase structure by expanding in powers of $ either 
the EoS or the effective potential and investigate its behaviour near the transition. 
In ref. [ E3 we chose the first possibility and determined the line of continuous 
transitions by requiring the vanishing of the square bracket of J2J at <& = 0. The 
pion mass was fixed at M = and only the fermion tadpole was expanded. The 
line of the 2 order phase transitions comes from the condition 



2 2 i ffV Ar , ( A g 2 N c 

The location of the TCP is determined by requiring in addition to m^ = also the 
vanishing of the coefficient of $ 2 in the square bracket of (J2J 
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Here ln(ci/2) = 1 — ^e + V an d V = Ih(Mqb/Mqf) gives a relation between M b 
and Mqf, the renormalisation scales of the bosonic and fermionic tadpoles. The 
parameters of the model are fixed at T = fi = 0. Choosing m q o = mjv/3 sw 312.67 
and requiring <&o = /w/2 for the solution of the EoS in the case of two flavors 
(Nf = 2) we get g = 6.72. A = 400 is fixed by requiring a good agreement between 

the location of the complex pole of the sigma propagator on the 2 nd Riemann sheet 
and the experimentally favoured mass and width of the sigma particle. Then one 
can solve the above equations and find the line of 2 order phase transitions, the 
first spinodal line and the location of TCP. This is shown in the l.h.s. of Fig. ^ 
The other lines in the figure are obtained numerically. In the r.h.s. we see the 
variation of the TCP and critical temperature with r/ and Mob- The study of the 
pole structure of the a propagator restricts the allowed range of Mqb to the right 
from the arrow, appearing on the horizontal axis of the figure. In this region no 
low scale imaginary poles appear besides the well known large scale tachyon, whose 
scale increases with decreasing values of the coupling constant A. In this region one 
can arrange very easily T c (/j, = 0) to occur in the region 150-170 MeV but the TCP 
stays robustly below 70 MeV. Previous effective model studies [|7j, and a recent 
renormalisation group improved investigation [ 111) , give similar values which are 
very low when compared with the CEP temperature of the QCD obtained in [|3] 
for physical number of quark flavors and masses. This difference could mean that 
the phase transition is not driven by the light degrees of freedom but by higher 
excited hadronic states. This possibility was supported by [ ^] when trying to 
reproduce lattice results on quark number susceptibilities and pressure with a non- 
interacting hadron resonance gas model. It was shown in [ I13| using renormalisation 
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group techniques that the inclusion of gluon degrees of freedom is important around 
T = 150 MeV. On the other hand we should keep in mind that we worked in the 
limit of an infinitely heavy strange quark, and that lattice studies revealed a very 
strong dependence on the values of quark masses [I14|. 

We now discuss a second method, in which we keep the ^-dependence of the 
mass gap M. The effective potential is obtained by integrating the functional form 
of EoS away from equilibrium. This analysis reveals an interesting property of 
the O(N) models in the large N limit: the vanishing of the effective coupling at the 
critical point [E3- To see this let us first note that <& = becomes a minimum of the 
potential as T — ► T c and so M — » 0. The fourth derivative of the effective potential 
at this point is Vgg($ = 0) = 4ff'(0) where H'(<5> 2 ) = ^(l- ^=^I F (M, m q )) - 



1g 2 N c K, T 2 lF(M,m q ) 
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Since the last term vanishes at T c one has to study the behaviour 
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Using a high temperature expansion Tb(M ) ~ j^ — ^~ 
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This means that along the line of 2 nd order phase transitions V"a vanishes. 
Still one should be able to provide a link between the effective potential of the 
model and Landau's theory of the tricritical point in which the location of the TCP 
is obtained from the condition of a vanishing quartic coupling. In order to achieve 
this, we first perform a high temperature expansion in 7J(<& 2 ) given by @ which 
yields 

H(^)=m 2 T + ^ + ^-2UM + -^ln C -f, (8) 
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where U = ££- an d K 

487T 
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In the 



following discussion for the sake of simplicity we omit in QSJ the one-loop fermion 
contribution Ip(M, m q ). We then rewrite the resulting iJ(4> 2 ) = M 2 relation as an 
equation which determines the <I> 2 dependence of M 



WM 2 + 2UM 



-4> 2 + K $4 
6 



0. 



(9) 



where W = 1 — ^^ In ^- and In J 2 - = i — 7. The solution of © approximates 
i/(<I> 2 ), so after substituting it into J2J it can be integrated to obtain an approxi- 
mation to the effective potential. Subsequent expansion in powers of $ 2 yields 
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A few remarks are in order here. First of all whether this potential has the right 

shape, and whether this approximation has the correct physical meaning depends 

on the actual values of the parameters. One can easily verify that for the parameters 

used in drawing Fig. U) W < in the temperature range of the 2 nd order phase 

transition. For increasing values of fi/T, k changes sign at 1.91 from negative to 

positive values. Since around TCP n has to be positive, this sign change restricts 

the location of the TCP to large values of [i: it must be in the region [i > 2T. 

Second, in the broken symmetry phase m 2 ^ < and it vanishes at T c meaning that 

the coefficient of the quartic term in <& vanishes also. Around T c from © one gets 

2 
m 2 ^ ~ T — T c . Therefore the scaling of the effective mass m c s = ^ryw ~ \T — T c \ 

gives the correct leading order critical exponent v = 1 for the O(N) model at large 

N . Third, to the expression in the square bracket in (|1L)|) the usual Landau type 

analysis applies. The location of the TCP is determined by the same conditions as 

in the analysis based on EoS, namely JfJJl and Q. To obtain the scaling exponent 

of the order parameter on the second order line we need only the first two terms 

in G3U- The minimum condition gives $ 2 = -^ =*> $ ~ (T - T c f , (3 = 1/2. 

At TCP one can set A = in (|1(J|1 and keep the sixth order term. The minimum 

2 
condition gives <I> 4 = —3^ =>■ <& ~ (T — Tc) 13 , f3 = 1/4. We verified numerically 

using the exact EoS that the value of (3 on the second order line and at the TCP 
agrees within the numerical accuracy with their mean field estimates given above. 
Refinement of the approximation is needed if one wants to enter the metastable 
region between the spinodals where m\ > and A < 0, but we think that never- 
theless the presented analysis captures essential features of the phenomenon. 



4. Conclusions 



In the 1/\A/V order of the large N approximation to the solution of the chiral quark- 
meson model we presented in the chiral limit an analytical determination of the 



T — jj, phase diagram of the chiral quark model 



phase diagram and of the location of the tricritical point. We have shown that with 
a reasonable approximation which works fine at high temperature, the effective 
potential of the theory can be related to mean field Landau type potential, from 
which the critical exponents can be determined. 

In order to assess to what extent a description of the phase diagram of QCD is 
possible using a low-energy effective model, one has both to improve our method of 
solving the model and to get as close as possible to the physical case of 2 + 1 flavors. 
The first issue necessitates taking into account the full momentum dependence of the 
self-consistent pion propagator, the implementation of a new resummation method 
and solving the Dyson-Schwinger equation for the fermions. The latter problem 
requires the extension of our method to the SU(3)l x SU(3)r case. 
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